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Aiming  at  the  general  problem  that  both  for  the  production  of  foundry  coke  made  by  anthracite  pulver¬ 
ized  coal  and  furnace  coke  made  by  tamped  coke  in  the  coke  oven,  namely  the  shrinkage  problem  for 
large-size  briquette  and  tamped  coke  in  the  pyrolysis  process,  a  mass  and  volume  integrated  measuring 
device  was  used  to  study  the  shrinkage  process  of  large-sized  briquettes  (volume  equals  103  cm3)  under 
constant  temperature  (450  °C,  550  °C,  650  °C  and  850  °C).  The  briquette  shrinkage  kinetic  model  was  also 
established  according  to  the  regulation  of  mass  loss  and  volume  shrinkage,  where  the  volume  shrink  cal¬ 
culation  method  of  tamped  coal  cakes  (13,280  mm  x  500  mm  x  4300  mm)  was  also  obtained.  Three  con¬ 
clusions  were  drawn  in  this  study:  1.  In  the  process  of  the  carbonization,  the  dehydration  and  degassing 
reaction  before  450  °C  as  well  as  polycondensation  reaction  at  650  °C  were  the  main  reason  causing  the 
briquette  shrinkage,  the  decomposition  reaction  at  550  °C  came  second,  while  the  hydrogen  evolution 
reaction  after  850  °C  has  the  least  effect  on  briquette  volume  shrinkage.  There  were  three  main  reasons 
causing  the  shrinkage  of  large-size  briquette:  the  decrease  of  coal  volume  caused  by  the  volatilization  of 
moisture  and  volatiles,  the  decrease  of  coal  volume  caused  by  the  pyrolysis  reaction  and  the  density 
increase  of  briquette  caused  by  the  decrease  of  coal  volume.  2.  The  briquette  shrinkage  kinetic  model 
was  represented  as  follows:  dVr/dt  =  V0[(l  -  x)d/J/df  -  A0e  £/RT(l  -  x)n/l];  There  was  a  positive  correlation 
between  the  volume  changing  parameter  /j  and  the  density  of  briquette,  and  the  volume  changing  param¬ 
eter  increased  with  the  reaction;  the  shrinkage  reaction  was  divided  into  three  stages;  with  the  reaction 
went  on,  the  reaction  order  and  reaction  rate  constant  decreased,  the  pre-exponential  factor  increased 
and  the  activation  energy  decreased  firstly  and  increased  then.  3.  The  shrinkage  of  the  tamped  coal  cakes 
from  the  large-scale  chamber  was  counted  in  the  calculated  unit  whose  thickness  was  100  mm.  The 
shrinkage  velocity  of  the  tamped  coal  cakes  in  carbonization  chamber  was  obtained  by  weighted  calcu¬ 
lating  the  shrinkage  parameters  of  each  volume  unit  at  different  temperatures.  The  shrink  regulations  of 
the  tamped  coal  cakes  in  the  carbonization  chambers,  which  were  in  agreement  with  industrial  produc¬ 
tion,  were  approximately  the  same  as  those  of  the  large-sized  briquettes.  However,  the  shrinkage  rate 
curve  of  the  tamped  coal  cake  in  the  carbonization  chamber  was  gentler  than  that  of  the  large-sized  bri¬ 
quette,  as  a  result  of  the  large  temperature  gradient  of  different  coal  in  carbonization  chamber. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

There  are  about  100  million  tons  of  high  rank  pulverized  coal 
burned  to  make  electricity  every  year,  making  resource  utilization 
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unreasonable  [1],  The  material  of  foundry  coke,  high-quality 
coking  coal  and  fat  coal,  is  turning  down.  As  a  result,  using  the  high 
rank  pulverized  coal  as  the  material  to  make  carbonization 
briquette  will  expand  the  type  of  coking  coal  and  relieve  the  envi¬ 
ronmental  pollution  at  the  same  time  [2],  Co-coking  is  commonly 
used  in  the  production  of  tamped  coke  in  the  large-scale  coke  oven, 
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Nomenclature 

V 

the  volume  of  briquette  (m3) 

i 

the  volume  unit  number  of  coal  material  in  the  carbon- 

W 

the  weight  of  briquette  (g) 

ization  chamber 

V 

the  volume  of  coal  in  carbonization  chamber  (m3) 

0 

the  physical  quantities  at  the  initial  time  point 

t 

carbonation  time  (s) 

t 

the  physical  quantities  at  the  time  t 

T 

carbonation  temperature  (K) 

f 

the  physical  quantities  at  the  end  of  the  reaction 

P 

the  volume  changing  parameter 

X 

the  conversion  rate  (%) 

Defined 

parameter 

k 

the  reaction  rate  constant 

Vt/Vo 

the  proportion  of  Vr  in  V0 

Ao 

the  pre-exponential  factor 

AVt/V0 

the  result  of  Vt/V0  at  high  temperature  minus  the  Vt/V0 

E 

the  activation  energy  (J/mol) 

at  its  lower  neighbor  temperature 

R 

ideal  gas  law  constant,  =8.314J/mol  K 

P 

the  proportion  of  the  shrinkage  degree  (A Vt/V0)  at  dif- 

n 

the  reaction  order 

ferent  temperature  zones  accounting  for  the  total  bri¬ 
quette  shrinkage  at  850  °C 

Subscripts 

a 

The  width  parts  of  tamped  coke  in  the  carbonization 
chamber  that  been  divided 

thus  the  shrinkage  characteristic  will  be  the  basis  to  guide  the 
production  of  large-scale  coke  oven  further.  Using  anthracite  or 
meagre  pulverized  coal  as  the  material  to  make  carbonization  bri¬ 
quette  and  tamped  coke  cleanly  is  an  important  technological 
approaches  to  use  high  rank  pulverized  coal  and  save  precious  cok¬ 
ing  coal  resource  scientifically  and  efficiently. 

The  shrinkage  characteristics  of  high-temperature  heat-treated 
briquettes  have  a  direct  effect  on  processes  such  as  the  discharge  of 
the  carbonized  briquette,  which  has  been  widely  studied  [3-5],  The 
properties  of  coal  and  carbonization  conditions  influence  briquette 
shrinkage.  The  properties  of  coal  include  the  coal  type  (e.g.,  volatile 
and  maceral),  binder  properties,  the  content  of  inert  material  and 
the  molding  density.  The  volume  shrinkage  rate  increases  linearly 
with  the  increase  of  coal  blending  volatiles  [6-8].  Of  the  three  pri¬ 
mary  types  of  macerals,  inertinite  promotes  briquette  shrinkage, 
vitrinite  has  no  obvious  relation  to  the  briquette  shrinkage,  while 
heating  vitrinite  promotes  briquette  expansion  [9,10];  binder  types 
have  a  significant  effect  on  the  volume  change  rate  of  the  bri¬ 
quettes,  and  the  high  ash  content  binder  increases  the  shrinkage 
level  of  the  first  shrinkage  peak  after  it  is  heated  and  solidified 
[9],  Adding  inert  materials,  such  as  coke  breeze  also  reduce  the  bri¬ 
quette  shrinkage  rate/contraction  coefficient  [11-13],  In  addition, 
the  briquette  shrinkage  rate  increases  with  the  compaction  density 
[14-16].  The  carbonization  conditions  include  the  heating  rate,  fur¬ 
nace  temperature  and  carbonization  time,  among  others.  The 
radial  shrinkage  rate  increases  with  the  heating  rate,  which  has 
little  effect  on  briquette  shrinkage  in  the  low  temperature  zone 
(approximately  500-750  °C),  while  the  effect  in  the  high  tempera¬ 
ture  zone  (approximately  750-960  °C)  is  significant  [17],  The  fur¬ 
nace  temperature  and  carbonization  temperature  influence  the 
briquette  shrinkage  primarily  by  heating,  which  has  an  effect  on 
fissure  formation  and  coke  density  [11], 

There  are  two  shrinkage  rate  peaks  during  the  coal  carboniza¬ 
tion  process.  First,  the  coal  undergoes  a  drastic  decomposition 
reaction,  precipitating  mostly  tar  and  heavy  hydrocarbon  gas;  the 
second  decomposition  process  occurs  in  the  solid  phase  and  pri¬ 
marily,  precipitates  hydrogen  [18],  Recent  studies  suggest  that 
the  free  radical  concentration  of  anthracite  increases  in  the  heating 
process,  and  the  free  radicals  combine  or  undergo  the  polyconden¬ 
sation  reaction  with  the  binder  [19].  The  shrinkage  is  caused  by  the 
free  radical  fragment  orientation  and  the  crystalline  ordered  array 
[17  .  The  coal  shrinkage  process  can  be  divided  into  three  stages: 
slow  contraction  stage,  transitional  contraction  stage,  and  rapid 
contraction  stage,  whose  activation  energy,  pre  exponential  factor 
and  rate  constant  increase  in  stages. 


The  shrinkage  of  coal  during  the  process  of  heat  treatment  affects 
its  pore  structure  [18],  reactivity  [19],  crack  formation  [20,21], 
strength  [22]  and  size  distribution  [23],  which  has  a  direct  relation¬ 
ship  with  gasification  characteristics  of  carbonized  briquettes. 
Recently,  the  research  of  powdered  coal  used  to  make  carbonized 
briquettes  for  the  gasification  industry  has  increased.  Currently, 
researchers  have  focused  on  single  coal  particles  or  small-sized  bri¬ 
quettes  (volume  smaller  than  103  mm3),  and  there  are  no  published 
studies  on  the  shrinkage  characteristics,  mechanism  and  kinetics  of 
large-sized  briquettes.  However,  studying  the  shrinkage  character¬ 
istics  of  large-sized  briquettes  is  imperative  for  the  industrial  use  of 
briquettes.  Thus,  a  mass  and  volume  integrated  measuring  device 
was  used  to  study  the  shrinkage  process  of  large-sized  briquette 
(volume  equals  103cm3)  under  constant  temperature(450  °C, 
550  °C,  650  °C  and  850  °C).  The  briquette  shrinkage  kinetic  model 
was  also  established  according  to  the  regulation  of  mass  loss  and 
volume  shrinkage,  where  the  volume  shrinkage  calculation  method 
of  tamped  coal  cakes  (13,280  mm  x  500  mm  x  4300  mm)  was  also 
obtained. 

2.  Experimental 

2.1.  Raw  materials 

Jincheng  anthracite  (JA)  and  Modified  coal-based  binder  (MB) 
was  used  in  this  experiment.  The  proximate  and  ultimate  analysis 
dates  are  shown  in  Table  1. 

2.2.  Mass  loss  measurement  of  large-sized  briquette  shrinkage 

The  anthracite  coal  and  binder  were  crushed  by  60  and  200 
orders  of  magnitude,  respectively.  The  anthracite  and  binder  were 
mixed  according  to  a  certain  proportion,  and  quantitative  amounts 
of  water  were  poured  into  the  mixture  and  stirred.  The  mixture 
was  put  into  a  cube  steel  mold  (100  mm  x  100  mm  x  100  mm) 
and  the  large  compound  briquettes  were  obtained. 

The  mass  and  volume  integrated  measuring  device  was  used  in 
this  experiment  during  the  pyrolysis  process  at  high  temperatures 
and  atmospheric  pressure  (Fig.  1 ).  The  measuring  device  contained 
as  quality  measurement,  volume  measurement  and  signal  process¬ 
ing  unit,  which  could  concurrently  obtain  the  measurements  of  the 
large-sized  briquette’s  mass  loss  and  volume  shrinkage;  the  sensi¬ 
tivity  of  mass  loss  and  volume  change  measurements  are  0.1  g  and 
0.1  mm  respectively  [24], 


Table  1 

Proximate  analysis  and  ultimate  analysis  of  raw  material. 
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Proximate  analysis  (wt.%) 

Ultimate  analysis  (wt.%) 

Mad 

Aad 

Vac, 

Q  af 

Hdaf 

Odaf 

Ndaf 

$daf 

Jincheng  anthracite  (JA) 

9.99 

11.44 

7.51 

92.74 

2.89 

2.49 

0.87 

1.01 

Modified  coal-based  binder  (MB) 

0.96 

13.46 

27.49 

83.49 

3.18 

12.38 

0.88 

0.07 

9  10  11  12  13  14  15  16 


Fig.  1.  The  mass  and  volume  integrated  measuring  device  during  the  pyrolysis 
process  at  high  temperature  and  ordinary  pressure  atmosphere,  1:  mass  sensor;  2: 
mass  measuring  rod;  3:  the  down  liquid  seal  unit;  4:  tray;  5:  retort;  6:  briquette;  7: 
furnace  body;  8;  Silicon  carbide;  9;  the  upper  liquid  seal  unit;  10:  air-inlet  duct;  11: 
temperature  sensor;  12:  exhaust  duct;  13:  volume  measuring  rod;  14:  the  light 
source;  15:  light  sensor;  16:  display  and  control  panel. 


Different  large-sized  briquettes  were  heated  to  450  °C,  550  °C, 
650  °C  and  850  °C  and  their  mass  loss  and  volume  change  was 
recorded.  The  temperature  was  kept  constant  and  the  experiment 
was  stopped  when  the  volume  change  and  mass  loss  were  less 
than  1  mm  and  0.1  g,  respectively.  The  regulation  of  the  volume 
change  and  mass  loss  of  the  briquettes  was  measured  at  different 
heating  rates. 


2.3.  Shrinkage  measurement  of  coal  in  carbonization  chamber 

The  industrial  experiment  for  the  shrinkage  model  was  per¬ 
formed  in  a  TJL4350D  stamping  coke  oven,  whose  size  parameters 
is  shown  in  Table  2.  Put  the  coal  in  the  coal  blending  ratio  that  as 
the  briquette’s  into  the  chambers.  The  coal  and  coke  positions  were 
measured  in  the  tamping  coke  oven  and  the  coal  volume  shrinkage 
rate  was  obtained. 


2.4.  The  evaluation  of  parameters 

2.4.1.  Sensitivity  analysis  [25] 

Sensitivity  analysis  is  commonly  employed  to  assess  that  in  the 
nonlinear  parameter  estimation,  and  the  nonlinear  parameter  esti¬ 
mation  is  usually  carried  out  by  using  the  least  squares  method  in 


Table  2 

Size  parameters  of  carbonization  chamber. 


Working 

Width  (pusher 

Width  (coke 

Width 

Height 

length  (mm) 

side,  mm) 

side,  mm) 

average 

(mm) 

(mm) 

13,280 

495 

505 

500 

4300 

order  to  find  the  global  minimum  of  the  following  objective 
function: 


(1) 


Sensitivity  analysis  is  applied  to  each  of  the  estimated  parame¬ 
ters  by  means  of  perturbations  of  the  parameter  value  (keeping  the 
other  parameters  in  their  estimated  values).  Perturbations  are 
preferably  done  in  the  range  of  ±20%.  For  each  perturbation  in 
the  parameter  values  the  objective  function  (the  volume  shrinkage 
rate  equation)  is  reevaluated  and  then  for  each  parameter  the  per¬ 
turbation  percentage  is  plotted  against  the  corresponding  value  of 
the  objective  function.  If  all  perturbations  in  all  the  parameters 
give  the  minimum  of  the  objective  function  with  their  original  val¬ 
ues  (0%  perturbation),  then  the  global  minimum  has  been 
achieved.  On  the  contrary,  if  at  least  one  parameter  does  not  give 
the  same  minimum  than  the  others  at  0%  perturbation,  that  means 
poor  nonlinear  parameter  estimation.  To  correct  this,  the  values  of 
the  wrong  estimated  parameters  are  re-determined  by  examining 
the  sensitivity  plot,  and  finally,  parameter  sensitivity  is  carried 
out  again  on  these  parameters  and  now  the  global  minimum  is 
guaranteed. 


2.4.2.  Residual  analysis 

Calculate  the  briquette  shrinkage  rate  difference  between  the 
experimental  value  and  optimized  value  to  investigation  the  prox¬ 
imity  of  them. 


Residual  = 


dV 

dt  exp 


dV 

dt  opt 


(2) 


2.4.3.  Parity  plots  analysis 

A  parity  plot  is  a  scatterplot  that  compares  expected  data 
against  experimental  data.  Each  point  has  coordinates  (x,  y),  where 
x  is  the  experimental  value  of  the  volume  shrinkage  rate,  and  y  is 
the  corresponding  optimized  value.  A  line  of  the  equation  y  =  x  is 
added  as  a  reference. 

3.  Results  and  discussion 

3.1.  Briquette  shrinkage  kinetic  analysis  during  the  carbonization 
process 

3.1.1.  Shrinkage  process  analysis  of  large-sized  briquettes 

The  briquette  shrinkage  curves  at  different  carbonization  tem¬ 
peratures  are  shown  in  Fig.  2a.  The  proportion  of  Vt  in  V0  is  repre¬ 
sented  by  VJVo.  The  Vf/V0  at  high  temperature  minus  the  VcjV0  at  its 
higher  neighbor  temperature,  then  the  A Vt/V0  is  obtained.  The  pro¬ 
portion  of  the  shrinkage  degree  ( AVJVq )  at  different  temperature 
zones  accounting  for  the  total  briquette  shrinkage  at  850  °C  was  p. 

As  can  be  seen  from  Table  3,  the  value  of  p  at  450  °C  is  large, 
dehydration  and  degassing  reaction  is  the  main  reaction  before 
450  °C  [24,26].  The  free  water  of  briquette  is  evaporated  by  heat, 
which  will  not  cause  the  volume  shrinkage  of  briquette  directly. 
Only  a  certain  volatilization  extent  of  free  water  will  lead  to  mea¬ 
surable  volume  shrinkage  [27].  The  decrease  of  coal  volume  caused 
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(C)  (d) 


Fig.  2.  Changing  curves  of  briquette  volume  with  time  (a)  shrinkage  curves;  (b)  shrinkage  rate  curve  (0-15  min);  (c)  shrinkage  rate  curve  (15-40  min);  (d)  shrinkage  rate 
curve  (40-140  min). 


Table  3 

Parameters  at  different  carbonation  temperatures. 


TV  C) 

(vyv0)  (%) 

(Avt/v0)  (%) 

P(%) 

450 

87.96 

12.04 

51.63 

550 

83.64 

4.32 

18.52 

650 

76.86 

6.78 

29.07 

850 

76.68 

0.18 

0.77 

by  the  volatilization  of  moisture  and  volatiles  and  the  decrease  of 
coal  volume  caused  by  the  pyrolysis  reaction  are  the  main  reason 
causing  the  shrinkage  of  large-size  briquette.  Decomposition  and 
condensation  reactions  carried  out  in  the  range  of  450-550  °C 
and  550-650  °C  separately.  The  decrease  of  coal  volume  caused 
by  the  pyrolysis  reaction  and  the  density  increase  of  briquette 
caused  by  the  decrease  of  coal  volume  are  the  main  reasons  caus¬ 
ing  the  shrinkage  of  large-size  briquette.  The  p  at  850  °C  was  0.77%; 
the  gas  precipitation  decreased  during  this  period,  which  was  the 
least  important  factor  causing  the  shrinkage. 

Fig.  2b-d  shows  the  briquette  shrinkage  rate  curves  under  dif¬ 
ferent  carbonization  temperatures.  The  instantaneous  shrinkage 
rates  of  the  large-sized  briquettes  were  infinite,  and  the  initial 
shrinkage  rate  was  calculated  according  to  the  reaction  values. 
When  the  reaction  was  stable,  the  shrinkage  rate  decreased  rapidly 
in  a  short  time  period.  A  period  later,  the  briquettes  were  in  the 
slow  shrinkage  stage,  and  the  shrinkage  rates  gradually  decreased. 

3.1.2.  Establishing  the  shrinkage  model 

Under  constant  temperature,  the  volume  of  the  briquette  made 
under  same  conditions  at  time  t(Vc)  was  only  related  to  the  initial 
volume  V0  and  the  conversion  x.  Under  the  same  initial  volume,  the 
larger  the  conversion,  the  smaller  the  Vt.  Therefore,  the  briquette 
volume  at  time  t(Vc )  had  a  positive  correlation  with  the  residual 
rate  1  -  x,  i.e„  Vtoc  Vo(l  —  x),  and  the  parameter  fl  was  introduced 
to  describe  the  volume  during  briquette  pyrolysis  process: 


Vt  =  V0{\-x)p  (3) 

where  /S  is  the  volume  changing  parameter  during  the  briquette 
pyrolysis  process,  which  varies  with  time  and  temperature. 

Under  constant  temperature,  the  relationship  of  the  volume 


change  rate,  conversion  and  volume  change  coefficient  is  obtained 
by  derivation  of  Eq.  (3): 

dVt/dt  =  V0[-(dx/dt)/J  +  (1  -  x)d/?/df] 

(4) 

where  combined  with  the  Arrhenius  equation 

k  =  A0e~E/RT 

(5) 

and  the  shrinking  core  model 

^  =  /<(!- x)" 

(6) 

the  shrinkage  kinetic  model  is  obtained 

dl/f/dt  =  V0[(l  -  x)dp/dt  -  A0e-E/RT(\  -  x)nf] 

(7) 

The  term  volume  changing  parameter  is  obtained  by  Eq.  (3),  and 
described  as  follows: 

P=Vt/[V  o(l-x)] 

(8) 

3.1.3.  Verification  of  the  shrinkage  model 

The  variation  in  the  regulation  of  the  briquette  volume  and 
mass  under  constant  temperature  was  obtained  by  experiment. 
The  theoretical  value  of  dVf/dt  was  obtained  through  Eq.  (7),  and 
by  derivation  of  the  equation  between  the  briquette  volume  and 
time,  the  experiment  value  of  dVJdt  was  obtained.  The  theoretical 
and  experiment  value  are  explained  in  Sections  3. 1.3.1  (calculation 
of  theoretical  value)  and  3. 1.3.2  (calculation  of  experimental  value) 
below  respectively. 
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Fig.  3.  Relation  graphs  of  ln(l  -  x)  and  ln(dx/dt)  under  550  °C:  (a)  the  diagram  of  whole  time  and  (b-d)  the  diagram  of  each  time. 


3. 1.3.1.  The  calculation  of  theoretical  value.  The  linear  equations  are 
obtained  by  taking  the  logarithm  on  both  sides  of  the  shrinking 
core  model  (Eq.  (4)). 

ln(dx/dt)  =  nln(l  -x)  +  lnk  (9) 

where  x  is  the  briquette  conversion,  x  =  (W0  -  Wt)/(W0  -  Wfj  (%), 
W0  is  briquette  mass  at  the  initial  time  point  (g),  Wt  is  briquette 
mass  at  the  time  t  (g),  Wf  is  briquette  mass  at  the  end  of  the  reaction 
(g),  k  is  obtained  by  the  Arrhenius  equation  (k  =  A0e~EIRT),  where  A0 
is  the  pre-exponential  factor  (s~]),  E  is  the  activation  energy  (J/mol); 
R  is  the  ideal  gas  law  constant  (J/(mol  K)),  T  is  reaction  temperature 
(K),  t  is  reaction  time  (s),  and  n  is  reaction  order. 

Eq.  (9)  was  used  to  draw  curves  of  ln(l  -  x)  on  ln(dx/dt),  using 
the  line  slope  and  intercept  to  obtain  the  reaction  order  n  and  coef¬ 
ficient  k  [24],  The  ln(l  -  x)  to  ln(dx/dt)  curve  at  550  °C  during  the 
briquette  pyrolysis  is  shown  in  Fig.  3a;  their  relationship  was  non¬ 
linear,  and  the  curve  was  divided  into  three  straight  lines  according 


to  the  approximate  reaction  stage.  The  graphs  can  be  seen  in 
Fig.  3b-d. 

The  related  graphs  of  ln(l  -  x)  and  ln(dx/dt)  at  other  tempera¬ 
tures  were  similar  to  the  graphs  at  550  °C.  The  reaction  order  n 
and  the  coefficient  k  at  different  temperature  stages  are  shown  in 
Table  4.  The  reaction  order  decreased  gradually  with  the  reaction 
at  the  same  temperature. 

The  reaction  order  decreased  with  the  reaction  time  at  the  same 
temperature,  and  that  is  similar  to  the  change  regulation  of 
briquette  pyrolysis  parameters.  Pitt  1 28]  figured  out  that  the  devol¬ 
atilization  of  coal  is  a  series  of  parallel  first  order  reaction.  A  heat¬ 
ing  process  is  needed  for  briquette  to  reach  the  reaction 
temperature  in  the  first  period,  and  there  are  kinds  of  parallel 
reactions,  thus  the  apparent  reaction  order  is  large.  With  the  time 
goes  by,  the  apparent  reaction  order  decreases  with  the  decrease  of 
reaction  kinds. 

The  reaction  rate  constant  decreases  with  the  reaction.  The 
reaction  rate  constant  reflects  the  shrinkage  speed  of  briquette  to 


Table  4 

Shrinkage  kinetic  parameters  of  briquette  at  different  reaction  temperature. 


Time  (h) 

r,  (°c) 

Hi 

o 

X 

R2 

r2  ("O 

n2 

k 2  X  105 

R2 

0-2 

2.580 

77.962 

0.947 

2.627 

83.950 

0.950 

2-5 

450 

0.625 

5.837 

0.952 

460 

0.788 

8.317 

0.987 

5-5.8 

0.165 

1.675 

0.985 

0.153 

1.699 

0.980 

0-1.25 

2.819 

139.634 

0.982 

2.877 

142.500 

0.983 

1.25-3 

550 

0.946 

24.750 

0.969 

560 

0.948 

24.777 

0.969 

3-4.5 

0.129 

3.766 

0.867 

0.132 

3.796 

0.872 

0-0.5 

2.838 

86.205 

0.981 

2.835 

86.403 

0.983 

0.5-2 

650 

1.306 

16.571 

0.992 

660 

1.306 

16.576 

0.992 

2-3 

0.694 

5.617 

0.995 

0.706 

5.760 

0.994 

0-0.5 

2.970 

419.736 

0.936 

2.997 

424.124 

0.938 

0.5-0.85 

850 

1.010 

42.391 

0.995 

860 

1.015 

42.514 

0.995 

0.85-1.4 

0.457 

14.373 

0.988 

0.462 

14.468 

0.988 
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Table  5 

Kinetic  parameters  of  briquette  shrinkage  at  different  reaction  temperatures. 


T  (°C) 

Time  (h) 

A0  x  10-3 

E  (J  mol  !) 

a 

0-2 

0.910 

56,077 

1 .248  -  0.0401nt 

450 

2-5 

5.694 

2618 

1.961  -0.1951nt 

5-5.8 

21.388 

6271 

0-1.25 

0.132 

11,585 

1.301  -  0.051  lnt 

550 

1.25-3 

3.687 

622 

2.054 -0.2211nt 

3-4.5 

13.636 

4524 

0-0.5 

0.936 

10,643 

1.299  -0.0591nt 

650 

0.5-2 

5.868 

216 

1.819 -0.1781nt 

2-3 

7.695 

5005 

0-0.5 

0.081 

9140 

1.313  -0.0631nt 

850 

0.5-0.85 

1.748 

2546 

1.707  -0.1941nt 

0.85-1.4 

3.518 

5800 

Table  6 

The  expression  of  briquette  volume  and  briquette  volume  shrinkage  rate  at  different 
reaction  temperatures. 


T  (°C) 

Vt/Vo 

(dV/dt)  x  103 

(m3s  ’) 

Experimental  value 

Theoretical  value 

450 

-0.0401nt  +  1.248 

0.040 

1 

t 

23.83t  +  3059 

550 

-0.0511nt  +  1.301 

0.051 

1 

t 

19.82t  +  2288 

650 

-0.0591nt  +  1.299 

0.059 

1 

t 

15.60t  + 1840 

850 

-0.063lnt  +  1.313 

0.063 

1 

t 

1 2.1 3t  +  1271 

some  extent.  At  the  initial  stage  of  the  reaction,  the  volatilization  of 
moisture  and  volatile  leads  to  rapid  volume  shrinkage;  the  reaction 
rate  constant  decreased  gradually  along  with  the  reaction  then. 
The  results  can  be  also  explained  by  the  content  in  Section  3.1.1. 
The  activation  energy  and  the  pre-exponential  factor  at  two  similar 
temperatures  were  approximately  the  same,  which  was  obtained 
by  the  Arrhenius  equation  through  the  calculation  of  the  k  value. 
Two  similar  temperatures  represented  by  Tj  and  T2  (T2  =  Tj  +  10) 


showed  a  coefficient  k  at  temperature  Tj  and  T2;  the  first  stages 
are  kj  a  and  k2  as  seen  below: 

ki=A0e-E/RT'  (10) 

k2=A0e-VRTi  (11) 

The  activation  energy  of  the  first  stage  at  temperature  T]  is  rep¬ 
resented  by  the  following: 

E  =  RTlT2/(T1-T2)ln(k1/k2)  (12) 


Bring  the  calculated  activation  energy  E  into  Eq.  (10);  the  corre¬ 
sponding  pre  exponential  factor  A0  is  obtained.  Bring  the  conver¬ 
sion  of  briquette  into  Eq.  (8),  the  volume  changing  parameter  p 
under  different  temperatures  are  obtained.  The  activation  energy 
E,  pre  exponential  factor  A0  and  the  volume  change  parameter  fl 
of  different  stages  and  temperatures  are  shown  in  Table  5. 

Under  the  same  temperature,  the  pre-exponential  factor 
increased  with  the  reaction.  And  that  is  similar  to  the  change  reg¬ 
ulation  of  briquette  pyrolysis  parameters.  From  the  collision  the¬ 
ory,  the  pre-exponential  factor  represents  the  collision  frequency 
of  the  molecules.  The  molecular  collision  frequency  increases  with 
increasing  temperature. 

The  activation  energy  of  large-size  briquette  decreased  firstly 
and  increased  then.  That  is  different  from  the  change  regulation 
of  briquette  pyrolysis  parameters  [29,30  .  Only  a  certain  volatiliza¬ 
tion  extent  of  free  water  will  lead  to  measurable  volume  shrinkage, 
thus  the  apparent  activation  energy  of  first  stage  is  large.  The 
shrinkage  of  briquette  is  mainly  caused  by  the  pyrolysis  of  coal 
in  the  second  and  third  stages;  as  a  result,  the  pyrolysis  activation 
energy  represents  the  shrinkage  activation  energy  in  some  extent. 
The  pyrolysis  activation  energy  increases  with  the  reaction,  and  so 
as  the  shrinkage  activation  energy. 

The  volume  changing  parameters  of  large-size  briquette 
increases  with  the  time.  And  there  was  a  positive  correlation 
between  the  volume  changing  parameter  /i  and  the  density  of  bri¬ 
quette  as  can  be  seen  from  Eq.  (13). 


V,  Vo  w  Vt  V,  V, 

V0(l  -x)*W0X  V0(l  -  x)  W0(l  -x)'XWt  Pt 


(13) 


0  20  40  60  80  100  120  140 


t/min 


Fig.  4.  Comparison  of  experimental  and  theoretical  values  at  different  temperatures. 
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The  pyrolysis  kinetic  parameters  of  briquette  was  calculated  on 
the  basis  of  mass  loss  change  regulation,  while  the  shrinkage 
kinetic  parameters  of  briquette  was  calculated  on  the  basis  of  the 
volume  change  regulation;  different  kinetic  model  basis  was  the 
main  reason  that  cause  the  difference  between  the  kinetic  param¬ 
eter  change  regulation  of  the  two  kinetic  model. 

Bring  the  shrinkage  kinetic  parameters  of  briquette  into  Eq.  (7) 
to  calculate  theoretical  value  of  briquette  shrinkage  rate,  and  the 
relational  expression  of  theoretical  value  was  obtained  (Table  6) 


3. 1.3.2.  The  calculation  of  experimental  value.  According  to  the 
briquette  shrinkage  curve  at  different  temperatures  in  Fig.  2a,  the 
briquette  volume  shrinkage  regulation  was  obtained.  The  experi¬ 
mental  value  of  dV/dt  is  acquired  through  the  differential  of  bri¬ 
quette  volume  shrinkage  regulation.  The  results  are  shown  in 
Table  6. 


The  experimental  and  theoretical  values  of  dVJdt  at  different 
temperatures  are  shown  in  Fig.  4,  and  the  results  are  consistent. 


3. 1.3.3.  The  evaluation  of  kinetic  parameters.  As  the  volume  chang¬ 
ing  parameter  (1  was  obtained  through  the  calculation  of  experi¬ 
mental  value  directly,  the  sensitive  analysis  was  only  practiced 
for  pre  exponential  factor  A0,  activation  energy  E  and  the  reaction 
order  n  at  450  °C,  550  °C,  650  °C  and  850  °C.  The  residual  analysis 
and  parity  plots  analysis  was  employed  to  compare  the  difference 
among  experimental  value  and  optimized  value  [25], 

As  the  kinetic  model  is  divided  into  three  sections,  the  sensitiv¬ 
ity  analysis  will  be  carried  out  in  three  sections  too.  For  each 
perturbation  in  each  parameter,  the  objective  function  was  evalu¬ 
ated  and  the  results  are  presented  in  Fig.  5.  The  mark  of  450  °C-1 
represents  the  first  stage  of  reaction  at  450  °C,  and  the  other  marks 
have  the  similar  meaning. 
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Fig.  5.  Sensitivity  analysis  of  kinetic  parameters  at  different  temperatures  and  during  different  time  period. 
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From  Fig.  5  it  is  clearly  seen  that  the  some  estimated  parame¬ 
ters  (such  as  the  curves  in  figure  450  °C-3)  are  the  optimum  since 
at  0%  perturbation  the  SSE  is  the  minimum,  which  is  the  condition 
to  assure  the  correct  values  of  parameters.  Some  parameters  in 
these  figures  (such  as  the  A0  curve  in  figure  450  °C-2)  do  not  show 
the  least  SSE  at  0%  perturbation,  and  then  optimized  procedures 
are  needed.  Fig.  6  finds  the  optimized  value  by  graphic 
visualization. 

The  optimized  value  of  kinetic  parameters  is  shown  in  Table  7. 

With  the  optimized  values  of  kinetic  parameters  sensitivity 
analysis  is  performed  again  and  the  results  are  shown  in  Fig.  7. 

The  residuals  analysis  and  parity  plot  analysis  of  four  tempera¬ 
tures  are  shown  in  Figs.  8  and  9  respectively.  The  shrinkage  kinetic 
model  of  large-size  briquette  in  the  pyrolysis  process  is  in  good 
agreement  with  the  experimental  value. 


3.2.  Industrial  amplification  of  the  shrinkage  model 

3.2.1.  Establishing  the  carbonization  chamber  shrinkage  model 

The  coal  material  from  the  carbonization  chamber  was  divided 
into  a  parts,  and  labeled  1,  2,  3, ...  i ....  a  from  left  to  right.  The 
initial  volume  of  coal  labeled  a  was  V'oi  (m3),  and  the  volume  of 
coal  at  time  t  was  V'ti  (m3).  In  industrial  production,  the  coal  in 
the  carbonization  chamber  is  heated  by  the  combustion  chamber 
on  its  two  sides,  which  was  similar  to  the  way  the  large-sized  bri¬ 
quettes  were  heated  in  this  study.  The  average  width  of  the  car¬ 
bonization  chamber  was  500  mm.  To  obtain  the  volume  units, 
the  carbonization  chamber  was  divided  into  10  parts  by  the  width, 
namely  the  parameters  were  adjusted  from  a  to  10;  thus,  the  coal 
labeled  i  and  a  —  i+  1  represented  the  volume  unit,  and  the  volume 
shrinkage  could  be  evaluated  by  the  shrinkage  of  the  briquettes 
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Fig.  6.  Secondary  sensitivity  analysis  of  kinetic  parameters  at  different  temperatures  and  during  different  time  period. 
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Table  7 

The  theoretical  value  and  optimized  value  of  briquette  shrinkage  parameters  at  different  reaction  temperatures. 


T("C) 

Time  (h) 

E  (J  mol 

-1) 

Ao 

X  to-3 

n 

Theoretical  value 

Optimized  value 

Theoretical  value 

Optimized  value 

Theoretical  value 

Optimized  value 

0-2 

56,077 

56,077 

0.910 

0.956 

2.580 

2.709 

450 

2-5 

2618 

2618 

5.694 

5.352 

0.625 

0.625 

5-5.8 

6271 

6271 

21.388 

21.388 

0.165 

0.165 

0-1.25 

11,585 

11,237 

0.132 

0.140 

2.819 

2.650 

550 

1.25-3 

622 

622 

3.687 

3.834 

0.946 

0.946 

3-4.5 

4524 

4524 

13.636 

13.636 

0.129 

0.129 

0-0.5 

10,643 

10,643 

0.936 

0.936 

2.838 

3.008 

650 

0.5-2 

216 

216 

5.868 

5.457 

1.306 

1.397 

2-3 

5005 

5005 

7.695 

6.849 

0.694 

0.756 

0-0.5 

9140 

10,237 

0.081 

0.072 

2.970 

3.297 

850 

0.5-0.85 

2546 

2546 

1.748 

1.748 

1.010 

1.010 

0.85-1.4 

5800 

5800 

3.518 

3.518 

0.457 

0.457 

%  Perturbation  %  Perturbation  %  Perturbation 


%  Perturbation 


%  Perturbation 


%  Perturbation 


Fig.  7.  Sensitivity  analysis  of  optimized  kinetic  parameters  at  different  temperatures  and  during  different  time  period. 


under  the  same  temperature  with  side-length  (100  mm).  The  bri¬ 
quette  shrinkage  rate  of  coal  in  the  carbonization  chamber  was 
obtained  by  weighting  the  shrinkage  parameters  of  volume  units 
under  different  temperatures. 

The  briquette  shrinkage  rate  could  be  obtained  from  Eq.  (4), 
where  the  volume  shrinkage  rate  of  coal  labeled  i  and  a  -  i  +  1  at 
time  t  is  shown  in  Eq.  (14): 

dV'./dt  =  Voi[(l  -  Xi)dft/d t  -  k,(  1  -  X,)"A]  (14) 


where  the  volume  changing  parameter  is 

Pi  =  V'a/[V'oi(\-Xi)]  (15) 

the  volume  shrinkage  rate  of  coal  in  the  combustion  chamber  is 

di/'t/dt  =  ^x^dv'./dr  (16) 

D  i 

and  the  volume  change  coefficient  of  coal  in  the  combustion  cham¬ 
ber  is 


10 
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Time/min  Time/min 


Fig.  8.  Residual  analysis  of  kinetic  parameters  at  different  temperatures. 


ft  =  vj 


v' 


E 

l 


(17) 


3.2.2.  Verification  of  the  carbonization  chamber  shrinkage  model 
According  to  the  volume  change  regulations  of  large-size  bri¬ 
quette  at  different  carbonization  temperatures  and  time  periods, 
the  volume  of  coal  labeled  i  and  a  -  i  +  1  at  different  temperatures 
at  different  carbonization  temperatures  and  time  periods  were 
obtained  and  the  estimated  value  of  volume  shrinkage  rate  for 
tamped  coke  in  champers  were  acquired  then.  According  to  the 


conversion  x,  the  Arrhenius  coefficient  k,  the  reaction  order  n,  vol¬ 
ume  changing  parameter  ft,  and  the  differential  of  volume  chang¬ 
ing  parameter  ^  of  large-size  briquette  at  different  carbonization 
temperatures  and  stages,  the  conversion  x ,,  the  Arrhenius  coeffi¬ 
cient  /q,  the  reaction  order  n„  volume  changing  parameter  /J,-,  and 
the  differential  of  volume  changing  parameter  ^  of  tamped  coke 
labeled  i  and  a  -  i  +  1  at  different  temperatures  were  calculated, 
and  the  theoretical  value  of  volume  shrinkage  rate  for  tamped  coke 
in  champers  were  acquired  further.  The  coal  and  coke  positions 
were  measured  in  the  tamped  coke  oven  and  the  volume  change 
regulations  of  tamped  coke  in  champers  was  calculated,  and  the 


Fig.  9.  Parity  plots  of  kinetic  parameters  at  different  temperatures. 
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position/mm 

Fig.  10.  Temperature  distribution  of  coal  in  carbonization  chamber. 


Table  8 

The  average  temperature  of  tamped  coke  labeled  i  and  a  -  i  +  1  during  different  time 
period. 


Time  (h) 

T,m 0  (°C) 

T2+9  (°C) 

T3+8  (°C) 

T4+7  (°C) 

T5+6  (°C) 

1 

550 

525 

175 

75 

50 

3 

800 

650 

450 

200 

75 

5 

825 

775 

600 

375 

200 

7 

885 

825 

700 

500 

300 

9 

950 

915 

775 

625 

450 

11 

990 

960 

860 

750 

600 

13 

1010 

1000 

940 

900 

820 

15 

1035 

1035 

1010 

990 

980 

Table  9 

VtIVo  of  briquette  at  different  temperatures  and  during  different  carbonation  time 
periods. 


V,/V o 


(°C) 

0-1  h 

1- 
3  h 

3- 
13  h 

13- 
15  h 

>15  h 

<350 

1 

1 

1 

_ 

- 

400 

1 

0.881 

0.881 

- 

- 

450 

0.01 07 12 

-  0.081 9t+ 1.0022 

0.881 

0.875 

- 

- 

500 

0.01 60t2 

-  0.1 045t+ 1.0024 

0.857 

0.848 

- 

- 

550 

0.0213I2 

-0.127U+ 1.0027 

0.833 

0.825 

- 

- 

600 

0.0457I2 

-  0.1 844t  + 0.9909 

0.805 

0.805 

- 

- 

650 

0.070212 

-  0.241 7t  + 0.9791 

0.777 

0.769 

- 

- 

700 

0.0882I2 

-  0.2727t  +  0.9827 

0.776 

0.768 

- 

- 

750 

0.106K2 

-0.3036t  + 0.9864 

0.774 

0.768 

0.768 

- 

800 

0.1 24U2 

-0.3346t  + 0.9900 

0.773 

0.767 

0.767 

0.767 

850 

0.168612 

-  0.3901 t  +  0.9965 

0.767 

0.767 

0.767 

0.767 

5*900 

0.243  It2 

-  0.451 6t+ 0.9753 

0.767 

0.767 

0.767 

0.767 

experimental  value  of  volume  shrinkage  rate  for  tamped  coke  in 
champers  are  acquired  further.  The  calculation  procedures  of 
established,  theoretical  and  experimental  value  were  explained 
in  Sections  3. 2. 2.1  (calculation  of  estimated  value),  3. 2. 2. 2  (calcula¬ 
tion  of  theoretical  value)  and  3.2.2. 3  (calculation  of  experimental 
value)  below  respectively. 

3.2.2. 1.  The  calculation  of  estimated  value.  Fig.  10  shows  the  temper¬ 
ature  distribution  of  coal  in  the  carbonization  chamber  after  differ¬ 
ent  carbonization  periods  [31].  The  average  carbonization 
temperature  of  coal  labeled  i  and  a  —  i  +  1  at  different  temperatures 
is  shown  in  Table  8. 

Table  9  shows  ratios  of  Vt/V0  for  large-size  briquette  at  different 
temperatures  and  time  periods  [32].  The  relational  expression 
between  the  shrinkage  volume  and  time  was  fitted  by  the  date 
within  1  h  at  different  temperatures.  When  the  temperature  was 
in  the  450-850  °C  range,  the  briquettes  shrank  completely  within 


Table  10 

The  Vi/Vo  of  tamped  coke  labeled  i  and  a  -  i  +  1  and  the  estimated  value  of  dV'/dt 
after  different  time  period. 


Time 

(h) 

V 1+10/ 
V'o 

V2  W 
V'o 

V'  3+8 / 

V  0 

vv7/ 

V'o 

V  5+6  / 

V  0 

V'  total  1 

V'0 

f  XU,- 

at  est 

(m3  min-1) 

0 

1 

1 

1 

1 

1 

1 

-4.75939 

1 

0.897 

0.906 

1 

1 

1 

0.961 

-2.37933 

3 

0.820 

0.881 

1 

1 

1 

0.940 

-3.09313 

5 

0.767 

0.769 

0.881 

1 

1 

0.883 

-2.69658 

7 

0.767 

0.768 

0.825 

1 

1 

0.872 

-2.30002 

9 

0.767 

0.767 

0.769 

0.825 

1 

0.825 

-2.97417 

11 

0.767 

0.767 

0.768 

0.806 

0.881 

0.797 

-2.30002 

13 

0.767 

0.767 

0.767 

0.768 

0.769 

0.767 

-1.18967 

15 

0.767 

0.767 

0.767 

0.767 

0.767 

0.767 

0 

100-140  min.  Thus,  the  volume  was  limited  by  the  carbonization 
time  in  the  range  of  0-3  h,  while  the  volume  was  the  ultimate  con¬ 
straint  under  corresponding  carbonization  temperature  in  the 
range  of  3-5  h.  Because  the  reaction  at  850  °C  has  less  effect  on 
the  briquette  shrinkage,  the  briquette  volume  at  temperatures 
higher  than  900  °C  was  similar  to  volumes  at  850  °C. 

The  following  are  examples  of  volume  calculation  procedures 
for  coal  in  the  carbonization  chamber. 

When  the  carbonization  time  was  1  h,  the  volume  of  tamped 
coke  labeled  i  and  a  -  i  +  1  are  calculated  directly  by  bring  the  car¬ 
bonization  time  into  relational  expressions. 

When  the  carbonization  time  was  3  h,  the  temperature  of  the 
carbonization  chamber  close  to  the  oven  wall  was  at  800  °C;  the 
temperature  volume  units  labeled  1  and  2  at  800  °C  were 
V'i+io  =  0.773  x  (0.2V'o),  and  the  average  temperature  volume  units 
labeled  2  and  9  at  450  °C  were  V2+a  =  0.881  x  (0.2\A0). 

The  volume  units  labeled  3  and  8  were  V'3+8  =  1  x  (O^V'o), 
the  volume  units  labeled  4  and  7  were  l/'4+7  =  1  x  (0.2Vo), 
and  the  volume  units  labeled  5  and  6  were  Va+6  =  1  x  (0.2V'o). 
Three  hours  later,  the  volume  of  the  briquette  in  the  carboniza¬ 
tion  chamber  was 

V  total =  (0.773  +  0.881  +  1  +  1  +  l)/5  =  0.931V'o _ 

The  V't/V'o  of  coal  labeled  i  and  a  -  i  +  1  after  different  time  peri¬ 
ods  were  calculated  in  the  same  way,  and  the  estimated  value  of 
volume  shrinkage  rate  ^esf  was  obtained  then,  the  results  are 
shown  in  Table  10. 

32.2.2.  The  calculation  of  theoretical  value.  According  to  the  bri¬ 
quette  volume  conversion  at  different  temperatures  [32],  the  con¬ 
version  Xj  of  coal  during  different  time  periods  and  at  different 
temperatures  is  shown  in  Table  11. 

The  examples  of  the  calculation  procedures  for  the  conversion 
of  coal  in  the  carbonization  chamber  can  be  seen  below. 

Fig.  10  shows  the  initial  temperature  of  the  oven  wall  liner  was 
1100  °C,  and  the  coal  was  placed  into  the  carbonization  chamber  at 
room  temperature.  Thus,  the  volume  unit  labeled  1  and  10  was  at 
550  °C  (1100  4-  2),  and  the  conversion  value  at  550  °C  was: 

x1+10  =  0.2295  x  ln3600  —  1.1135  =  0.766 

the  volume  unit  labeled  2  and  9  at  500  °C  had  a  conversion 

value  of 

x2+g  =  0.1954  x  ln3600  -  0.9237  =  0.676, 

the  volume  unit  labeled  3  and  8  at  250  °C  had  a  conversion  of 
x3+8  =  0.36, 

the  volume  unit  labeled  4  and  7  at  150  °C  had  a  conversion 
value  of  x4+7  =  0.15, 

and  the  volume  unit  labeled  5  and  6  at  50  °C  had  a  conversion  of 

x5+6  =  0. 


12 


Q.  Wang  et  al./Fuel  138  (2014)  1-14 


Table  11 

x,  at  different  temperatures  and  during  different  time  period. 


Temperature  (°C) 

X, 

Fitting  time  (h) 

0-1  h 

1-3  h 

>5  h 

50 

0 

0.10331nt- 0.5311 

0.13101nt  -  0.7303 

1.5 

100 

0.07 

0.1 1 1 71nt  —  0.5685 

0.13471nt  -  0.7534 

1.7 

150 

0.15 

0.12011nt- 0.6058 

0.14381nt  -  0.7765 

2.0 

200 

0.26 

0.12851nt- 0.6431 

0.15021nt- 0.7996 

2.5 

250 

0.36 

0.13691nt  -  0.6805 

0.15661nt- 0.8227 

3.0 

300 

-0.04331nt  +  0.4045 

0.14531m-  0.7178 

0.16301nt  -  0.8458 

4.0 

350 

0.02491nt  +  0.0250 

0.15361m-  0.7551 

0.16941nt  -  0.8689 

6.5 

400 

0.093  lint-  0.3545 

0.17041nt- 0.8300 

0.18221nt  -  0.9151 

6.2 

450 

0.16131nt- 0.7340 

0.18711nt- 0.9044 

0.19501nt- 0.9613 

6.0 

500 

0.19541nt  -  0.9237 

0.20391nt  -  0.9791 

0.20781nt-  1.0075 

4.5 

550 

0.22951nt-  1.1135 

0.22061nt 

1.0537 

3.0 

600 

0.18731nt- 0.8600 

0.19161m- 0.8801 

0.19941nt  -  0.9361 

3.8 

650 

0.14521nt- 0.6066 

0.16261nt  -  0.7065 

0.17811nt- 0.8185 

4.5 

700 

0.15661nt- 0.7225 

0.17051nt- 0.7065 

0.18211nt- 0.7905 

4.0 

750 

0.16791nt- 0.7110 

0.17831nt  -  0.7065 

0.18611nt- 0.7625 

3.5 

800 

0.17931nt- 0.6995 

0.186Unt- 0.7065 

0.19001nt  -  0.7345 

2.4 

850 

0. 19401nt-  0.7065 

1.8 

5*900 

0.19991nt- 0.6925 

1.5 

When  the  carbonization  time  was 

3  h,  the  calculation  proce- 

0.8  - 


0.6 


0.4 


0.2 


0.0  - 

2  4  6  8  10  12  14  16 

Time/h 

Fig.  11.  The  conversion  x,  of  tamped  coke  labeled  i  and  a  -  i  +  1  after  different  time 
period. 


dures  for  the  conversion  of  coal  were  the  following. 

Fig.  10  shows  the  temperature  of  the  oven  wall  was  approxi¬ 
mately  800  °C  after  3  h.  Thus,  the  volume  unit  labeled  1  +  10  was 
at  800  °C.  This  part  is  under  550  °C  before  1  h,  and  whose  conver¬ 
sion  value  is  0.7658,  which  equaled  the  reaction  time  t  =  0.9838  h 
(0.1 7931nt  -  0.6995  =  0.7658)  at  800  °C;  at  this  point,  the  coal 
had  completed  the  reaction,  and  the  conversion  at  this  stage  was 
zero.  The  conversion  of  coal  labeled  i  and  a  —  i  +  1  during  different 
time  period  was  calculated  in  the  same  way,  and  list  in  Fig.  11. 

As  the  same  calculation  method  for  the  conversion  of  the 
tamped  coke  labeled  i  and  a  -  i  +  1,  the  reaction  order  n„  the  Arrhe¬ 
nius  coefficient  ki}  volume  changing  parameter  and  the  differen¬ 
tial  of  volume  changing  parameter  ^  of  tamped  coke  labeled  i  and 
a  —  i  +  1  at  different  temperatures  were  calculated  according  to  the 


Time/h 


Time/h 


Fig.  12.  The  reaction  order  n,  (a),  Arrhenius  coefficient  k,  (b),  volume  changing  parameter  /?,  (c),  and  the  differential  of  volume  changing  parameter  ^  (d)  of  tamped  coke 
labeled  i  and  a  -  i  +  1  after  different  time  period. 
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Table  12 

The  dV'i/dt  of  tamped  coke  labeled  i  and  a  —  i  +  1  and  the  theoretical  value  of  dV',/dt  after  different  time  period. 


Time  (h) 

X  to-6  (m3-  min-1) 

i  =  1 

i  =  2 

i  =  3 

i  =  4 

i  =  5 

-j-  x  10  3  (m3  min  !) 

d  t  the 

1 

-1221.397 

-396.587 

0 

0 

0 

-1.618 

3 

-1031.013 

-1085.889 

-947.527 

0 

0 

-3.064 

5 

-420.343 

-896.133 

-1392.138 

0.017 

0 

-2.708 

7 

0.011 

-11.706 

-888.253 

-1393.28 

0.046 

-2.293 

9 

0.228 

-3.997 

-451.807 

-1542.094 

-948.212 

-2.946 

11 

0.171 

0.228 

-12.049 

-452.092 

-1835.151 

-2.299 

13 

0 

-0.114 

-3.997 

-305.449 

-900.644 

-1.210 

15 

0 

0 

0 

0.057 

-23.812 

-0.024 

Time/h 

(a) 


£ 

a 

© 

x  -« 


2.0 
1.5  - 
1.0 
0.5 
0.0 
0.5 
1.0 
1.5 
2.0 


□  theoretical  value 
o  estimated  value 


0  2  4  6  8  10  12  14  16 

Time/h 

(b) 


Fig.  13.  Comparison  of  the  theoretical,  estimated  and  industrial  value  for  the  shrinkage  rate  of  tamped  coke  in  carbonization  chamber:  (a)  shrinkage  rate  curves  and  (b) 
residual  analysis. 


kinetic  parameters  of  large-size  briquette  at  the  same  carboniza¬ 
tion  [32]  and  the  carbonization  temperatures  and  stages  of  the 
tamped  coke.  The  results  are  shown  in  Fig.  12. 

Put  the  kinetic  parameters  of  tamped  coke  labeled  i  and  a  -  i  +  1 
into  Eq.  (14),  the  theoretical  value  of  volume  shrinkage  rate  for 
tamped  coke  after  different  carbonization  time  was  obtained. 

When  the  temperature  was  lower  than  100  °C,  the  coal  had  not 
undergone  the  moisture  evaporation  reaction,  the  volume  of  coal 
did  not  change,  and  the  volume  shrinkage  rate  was  set  to  zero; 
when  the  temperature  was  higher  than  900  °C,  the  volume  of  the 
coal  change  little,  and  the  volume  shrinkage  rate  was  set  to  zero 
(see  Table  12). 

3.2.23.  The  calculation  of  experimental  value.  The  coal  and  coke 
positions  were  measured  in  the  tamped  coke  oven  and  the  volume 
change  regulations  of  tamped  coke  in  champers  was  calculated, 
and  the  experimental  value  of  volume  shrinkage  rate  for  tamped 
coke  in  champers  were  acquired  further  (Fig.  13a).  Residual  analy¬ 
sis  is  employed  to  compare  the  difference  among  the  estimated 
value,  theoretical  value  and  experimental  value. 

The  variation  regulation  of  the  theoretical,  estimated  and  exper¬ 
imental  value  for  the  carbonization  chamber  coal  shrinkage  rate  is 
shown  in  Fig.  13.  All  three  values  were  in  agreement.  The  shrink 
regulations  of  tamped  coal  cakes  in  carbonization  chambers  are 
consistent  with  the  large-sized  briquettes. 

However,  the  shrinkage  rate  curve  of  the  tamped  coal  cakes  in 
the  carbonization  chamber  was  smaller  than  for  the  large-sized 
briquettes  because  of  the  large  temperature  gradient  of  different 
coal  types  in  the  carbonization  chamber.  Specifically,  the  instanta¬ 
neous  shrinkage  rate  of  the  coal  material  at  high  temperatures  in 
the  carbonization  furnace  rose  by  6  x  1 0  3  m3/min  from  zero  at 
the  initial  time  point.  When  the  reaction  was  stable,  the  shrinkage 
rate  decreased  rapidly.  The  dehydration  and  degassing,  thermal 


decomposition,  and  the  condensation  reactions  occurred  during 
the  pyrolysis  process,  and  the  shrinkage  rate  increased  gradually; 
the  gas  precipitation  amount  decreased  in  the  later  reaction  stage, 
and  the  shrinkage  rate  gradually  returned  to  zero.  However,  the 
shrinkage  rate  curve  of  the  tamped  coal  cake  in  the  carbonization 
chamber  was  smaller  than  for  the  large-sized  briquettes  because 
of  the  large  temperature  gradient  between  the  different  coal  types 
in  the  carbonization  chamber. 

4.  Conclusion 

(1)  In  the  process  of  the  carbonization,  the  dehydration  and 
degassing  reactions  before  reaching  450  °C  and  the  polycon¬ 
densation  reaction  in  the  range  of  550-650  °C  were  the  pri¬ 
mary  causes  of  briquette  shrinkage  during  the  carbonization 
process;  the  leakage  of  gas  decreased  after  850  °C,  which  had 
the  least  effect  on  the  briquette  volume  shrinkage.  There 
were  three  main  reasons  causing  the  shrinkage  of  large-size 
briquette:  the  decrease  of  coal  volume  caused  by  the  volatil¬ 
ization  of  moisture  and  volatiles,  the  decrease  of  coal  volume 
caused  by  the  pyrolysis  reaction  and  the  density  increase  of 
briquette  caused  by  the  decrease  of  coal  volume. 

(2)  The  briquette  shrinkage  kinetic  model  was  represented  as 
follows:  dVf/dt  =  V0[(l  — x)d/3/dt  —  A0e  E/RT(1  — x)ny8];  There 
was  a  positive  correlation  between  the  volume  changing 
parameter  p  and  the  density  of  briquette,  and  the  volume 
changing  parameter  increased  with  the  reaction;  the  shrink¬ 
age  reaction  was  divided  into  three  stages;  with  the  reaction 
went  on,  the  reaction  order  and  reaction  rate  constant 
decreased,  the  pre-exponential  factor  increased  and  the  acti¬ 
vation  energy  decreased  firstly  and  increased  then.  The 
pyrolysis  kinetic  parameters  of  briquette  was  calculated  on 
the  basis  of  mass  loss  change  regulation,  while  the  shrinkage 
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kinetic  parameters  of  briquette  was  calculated  on  the  basis 
of  the  volume  change  regulation;  different  kinetic  model 
basis  was  the  main  reason  that  cause  the  difference  between 
the  kinetic  parameter  change  regulation  of  the  two  kinetic 
model. 

(3)  The  shrinkage  of  the  tamped  coal  cakes  from  the  large-scale 
chamber  was  counted  in  the  calculated  unit  whose  thickness 
was  100  mm.  The  shrinkage  velocity  of  the  tamped  coal 
cakes  in  carbonization  chamber  was  obtained  by  weighted 
calculating  the  shrinkage  parameters  of  each  volume  unit 
at  different  temperatures,  and  the  volume  shrinkage  calcula¬ 
tion  method  of  the  tamped  coal  cake 
(13,280  mm  x  500  mm  x  4300  mm)  was  obtained.  The 
shrink  regulations  of  the  tamped  coal  cake  in  the  carboniza¬ 
tion  chamber  were  consistent  with  the  large-sized  bri¬ 
quettes.  However,  the  shrinkage  rate  curve  of  the  tamped 
coal  cake  in  the  carbonization  chamber  was  smaller  than 
the  large-sized  briquette  because  of  the  large  temperature 
gradient  of  the  different  coal  types  in  the  carbonization 
chamber. 
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